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CHAPTER 1 


INTRODDCTIOH 

Computational Geometry , as it stands today , is con- 
cerned with the computational complexity of geometri? prob- 
lems within the framework of analysis of algorithms. A 
large number of applications areas such as pattern recogni- 
tion , computer graphics , image processing , operations 
research , statistics , computer aided design , robotics , 
etc. , have been the incubation bed of the discipline since 
they provide inherently geometric problems for which effi- 
cient algorithms have to be developed. These problems 
include the Euclidean travelling salesman , minimum spanning 
tree , linear programming , and hosts of others. Algo- 
rithmic studies of these and other problems have appeared in 
the scientific literature with an increasing intensity in 
the past two decades and a growing number of researchers 
have been attracted to this discipline , christened "Compu- 
tational Geometry" in a paper by Shamos [1] in 1976. 

According to the nature of the geometric objects 
involved , we can identify basically five categories into 
which the entire collection of geometric problems can be 
conveniently classified , i.e. , convexity , intersection , 
geometric searching , proximity , and optimisation. 

Since our research area is concentrated around the 



proximity problem , we shall briefly describe the proximity 
problem and the work done in this area in the next few sec- 
tions of this chapter, 

1. EBQXIMITY PROBLEMS 

Geometric objects , such as points and circles , are 
used to model physical entities in the real world. In some 
cases we would like to have access to a suitable neighbour- 
hood of the objects. For instance , in air traffic control 
we wish to keep track of the closest two aircrafts ; when 

aircrafts are modelled as points moving in space , we want 

to find the closest pair of points at a certain point in 

time. We shall first list a number of problems , some of 

which may appear unrelated , and describe a geometric con- 
struct , called a Voronoi diagratm , which can be used to 
solve these problems within the same order of time spent for 
computing the diagram. 

1.1. Basic Frckimlty prcbleias 

Problem 1 (Closest peir > > 

Given n points in the plane , find two points that are 
closest. 

It is obvious that the generalisation of the problem in 
k dimensions , k >= 1 , can be solved in 0(k n^) time by 
computing all interpoint distances. In one dimension , we 
can solve the problem easily in 0(n log n) time by a prelim- 
inary sorting. It turns out that sorting does not general- 
ise to higher dimensions. Using the divide-and-con<auer 



technique , Bentley and Shames [2] showed that 0(n log n) 
time is sufficient to solve this problem in dimensions k >= 
1 , and the time bound is optimal. An average case study of 
this problem is presented in [3] where an optimal average 
case algorithm is described. 

Problem 2 (All nearest neighbours! s 

Given n points in the plane , find for each point a 
nearest neighbour (other than itself). 

Problem 3 CEucl idean minimum spanning tree « EHSTls 

I 

Given n points in the plane , find a tree that inter- 
connects all the points with minimiam total edge length. 

This problem has an obvious application in computer 
networking where we want to interconnect all the computers 
at minimum cost. This formulation , however , forbids the 
addition of extra points. If additional points , called 
Steiner points , are allowed , the problem becomes the 
minimal Steiner tree problem , which has been shown to be 
NP-hard £4] . Note also that the EMST problem can be cast as 
a graph-theoretical problem , in which the weight of each 
edge is the distance between the two terminal vertices of 
the edge. In [5] several spanning tree algorithms have been 
illustrated. In general , the minimum spanning tree problem 

O 

for a graph with n vertices requires J\-(n ) time , for the 
minimum weight edge must be in the tree and there are O(n^) 
independent weights in the input ; however , the Euclidean 
metric properties can be exploited so as to solve the EMST 
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problem in 0(n log n) time. 

Problem 4 <Tr iangulation) s 

Given n points in the plane , construct a planar graph 
on the set of points such that each face within their convex 
hull is a triangle. 

This problem arises in numerical interpolation of 
bivariate data where the function values are known at irreg- 
ularly spaced points , and in the finite element method. A 
triangulation of these n points can be used to approximate 
the function value at a new point as the interpolation of 
the function values at the vertices of the triangle contain- 
ing the new point. 

Problem 5 (Nearest neighbour search! s 

Given n points in the plane , with preprocessing 
allowed , find the nearest neighbour of a querry point. 

This problem , also known as the "post office problem" > 
arises in pattern classification [6] where the nearest 
neighbour decision rule is used to classify a new sample 
into the class to which its nearest neighbour belongs to , 
and in information retrieval where the record that best 
matches the querry record is retrieved [73 . 

Problem 6 (k Nearest neighbours search) t 

The same as Problem 5 except that the k nearest neigh- 


bours are sought. 
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l.Z. IQE YPRQRQI DIAGRAM 

The above problems can be solved efficiently by using 
Voronoi diagram. Given a set S of n points {p^^ » ^2 ’ ' ‘ ’ ’ 
p^} , the Voronoi diagram [83 of S , denoted by Vor (S) , 

pariiitions the plane into n "equivalence" classes , each of 
which corresponds to a point. Specifically , the 
equivalence class corresponding to p^^ is the Voronoi polygon 
V (p^) , which is formally defined as V (p^) = {f i r in 
and d (r , Pj^) <= d (r , Pj) , j <> i}. In other words , V 
(p^) is the locus of points that are as close to p^^ as any 
other point of S and can also be defined as the intersection 
of the half planes (Pj^ , p^) where H (Pj^ , pj) is the 

half plane determined by the perpendicular bisector of p^^p^ 
and containing Pj^ . Thus , the Voronoi diagram of a set of 
n points is just a collection of n Voronoi (convex) polygons 
, one for each point. The diagram is also called Thiassmn 
polygons [93. The properties of Voronoi diagram have been 
discussed in details in [103 , [113. The straight-line dual 
of the Voronoi diagram is a triangulation of S. The tri- 
angulation is also known as Dalaunay triangulation and 
Dirichlet tessellation [123. The Voronoi diagram of a set 
S of n points in the plane can be constructed in 0(n log n) 
time , which is optimal. 

Once the Voronoi diagram ts available , the closest 
pair problem , the all-nearest-neighbour problem , and the 
triangulation problem can all be solved in Ofn^ 
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Obtain the straight line dual graph by scanning each edge of 
the Voronoi diagram ; since the dual graph is a triangula- 
tion and the total number of edges in Vor (S) is 0(n) , the 
process takes 0(n) time. The closest pair is identified 
with an edge of the triangulation and , similarly , the 
nearest neighbour of each point is given by an edge of the 
triangulation ; therefore , both problems can be solved in 
0(n) time. It has been shown [11] that the EMST is a sub- 
graph of the delaunay triangulation. So the EMST problem 
can also be solved in additional 0(n) time using the algo- 
rithm of Cheriton and Tar jan [6] . As for the nearest neigh- 
bour search problem , all we need to do is to find the Voro- 
noi polygon in which the new point lies. The search is 
therefore a point location problam [13] , [14] and can be 
carried out in O(log n) time. 

z. csamisLiim ssi im thesis 

In this thesis we have looked at various kinds of 
proximity problems. In the second chapter of this thesis , 
a new 0(n log n) algorithm is proposed for the closest-pair 
problem. This algorithm is based on a novel variation of 
the Sweep-line paradigm. Here we use a pair of lines whose 
distance from each other varies dynamically to sweep the 
point set. 

In the third chapter a type of generalisation of Voro- 
noi diagram called Kth degree Voronoi diagram is discussed. 
A new dynamising algorithm is given whereby a Kth degree 
Voronoi diagram can be updated efficiently. An easy-to- 
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implement incremental algorithm for constructing a Kth 
degree Voronoi diagram is also proposed in this chapter. 

In the fourth chapter we have dealt with the shortest- 
path problem inside a simple polygon containing a constant K 
number of obstacles. 

Finally , in chapter five we have concluded this thesis 
, mentioning some related open problems. 
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CHAPTER 2 

A NEfV ALGORITHM FOR THE CLOSEST PAIR PROBLEM 


1. INTRODPCTION 

The CLOSEST PAIR problem in compu-babional geometry is 
the following *• , 

Given a set of n points in the plane find a mutu- 
ally closest pair. 


Apart from its intrinsic interest , an efficient solution to 
this problem would be useful , for instance , in air-traffic 
control . 

Worst-case time optimal O(niogn) algorithms for this 
problem have been given baaed on (a) the Divide-and-Conquer 
approach , (b) the construction of the Voronoi diagram of 
the given point set [53- 

In this paper , we describe a new algorithm for the 

problem , based on a novel variation of the well known 

« 

sweep-line paradigm in computational geometry. 

The rest of the paper is organised as follows : In Sec- 
tion 2 we discuss a lower bound for the problem in the alge- 
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braic decision tree model. . Section 3 contains a brief 
description of the sweep-line paradigm , and a detailed dis- 
cussion of our variation of it , which we call the sweep- 
rectangle method. To motivate our solution to the problem , 
in Section 4 we discuss the case where all the points lie on 
a straight line. The general case is discussed in Section 6 
and in the next section we discuss the case of d-dimensions . 
We summarise in Section 7 and indicate future directions of 
research. 

2. LQHEK JXSm 

The ELEMENT UNIQUENESS problem is to decide if n given 
real numbers are all distinct. We can transform this prob- 
lem in linear time to the closest pair problem by consider- 
ing the set of real numbers {xl , x2 , ... , xn} as n points 
on the x-axis. Clearly , the elements are distinct if the 
distance between a closest pair is non-zero. Since in the 
algebraic decision tree model [1] [5] any algorithm that 
determines whether the elements of a set of n real numbers 
are distinct requires -fb (nlogn) tests , a lower bound of 

(nlogn) is established for the CLOSEST PAIR problem. 

2. SWEEP-BECTAHSIiE TECHNIQUE 

The sweep- line paradigm is best understood by means of 
an example to which it has been successfully applied. To 
report all intersecting pairs among a set of n line segments 
in the plane , we sweep the set from -(inf) in the 'left to 
+(inf) in the right by a straight line moving in a direction 




orthogonal to itself . The computation is supported hy two 
dynamic data structures : The sweep- line status and the 
event point schedule. The former keeps record of the set of 
line segments which currently intersect the sweep-line and 
the latter maintains a list of points called event points at 
which all computations (like testing for intersection , 
deletion & insertion of line segments etc.) are done [4] [5] 



In our variant of this technique , two parallel lines 1 

I 

and 1 , perpendicular to the x-axis , and at a distance d 

from each other which varies dynamically , is swept across 
the point sot from left to right as shown in the above fig- 
ure. (The two lines are the two sides of an infinitely long 
strip , and hence the name sweep- rectangle) . Instead of the 
sweep-line status , here we have the sweep-rectangle status 
, which maintains the points currently inside the sweep- 
rectangle. The event point schedule is a lexicographically 
sorted queue of the given points. The sweep- rectangle 
status is updated at the event points. The update involves 
insertion of a now event point in the sweep-rectangle status 



12 


and deletion of the points going outside it. This update 
over , distances are computed between the newly inserted 
point and its nearest neighbour candidates inside the 
sweep- rectangle. A shortest distance pair among these is 
used to check if a closest pair of the points seen so far 
needs updating. 

The width d of the rectangle decreases as it sweeps the 
point set. Initially it is chosen to be infinite and gets 
updated as the sweep goes on. At any moment of time , d 
gives the distance between a closest pair among the points 
seen so far. The elegance of this technique lies in the 
fact that the number of distance computations for each newly 
inserted point is constant. We can show that when a new 
point is inserted only the distances between the new point 
and a constant number of existing points are to be computed 
(discussed in details below). These constant number of 
points for each newly inserted f. point are called as the 
"nearest neighbour candidates" for that new point. 

Qasi dlmmalQiifll ssas& 

Let us first try to develop an algorithm for the one 
dimensional case , using the above technique. W.l.o.g. we 
may assume that all the points lie on the x-axis. In this 
case the sweep-rectangle degenerates to an interval on the 
x-axis. It sweeps the point set from -(inf) to +(inf). A 
formal algorithm for this case is given beloW. 
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Procedure SMeep-reci—one—dim; 

{It -finds a closest pair (p^ql of a set of H points 
lying on the x—axisl 
begin 

Sort the H points and place them in the 
queue Ef 

S s=> cp f /* S is the rectangle status »/ 
d s* (inf)f 

If the no. of points in E is more than one 
begin 

Extract the first tmo points p and q from E , find 
the distance between them and assign it to df 

5 •. r 5 u f k.cv I i 

»hile <E <> 4^ } do 

* 

begin 

Extract from E the first element pi 
Insert pi in Sf 

Delete from S all the points having 

abscissae less than pi — df 

Compute the distance of pi from the 

other points in Sf 

Find . the minimum of the above 

distances (say dllf 

Let ql be a point at a distance 

dl from pi f 

d min (dl , d)f 

If d = dl then /* update p & q »/ 
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begin 

P f= Pif 
q ql? 
end; 

end; 

end; 

Output d & (p f q} ; 

end; 


HaiA atrttgtu re a 

The data structure implementing the event point 
schedule has to support onl/ the operations MIN (E) , which 
determines the smallest element in E and deletes it. We can 
use a linear list implementation of E which supports the 
above operation in 0(1) time. 

The data structure for the sweep- rectangle status is 
determined by the following lemma: 

Lemmal: At any point of time , the rectangle contains at 

most two points. 

Proof: The proof of the above lemma follows from the 

observation that d is the distance between a 
closest pair of the points met so far. 

So , we can use a list having two elements to implement 
the sweep- rectangle status. This supports insert and delete 
operations in 0(1) time. 



Theoreml: The above algorithm correctlv finds a mutually 
closest pair for a one dimensional distribution of 
points . 

Proof: The proof follows from the fact that at any moment 

of time , d gives the distance between a closest 
pair and (p,q) gives the two points forming a 
closest pair for all the points to the left of the 
current event point. 

Theorea2: The above one dimensional closest pair algorithm 
takes OCnlogn) time and uses 0(n} space. 

Proof: There are 0(n) points in the point set and by 

Lemmal , for each newly inserted point , we have to 
make 0(1) distance computations. 0(n) insertions 
take O(nlogn) time and 0(n) deletions take 0(n) 
time. The initial sorting step takes O(nlogn) 
time. So , the overall time complexity is 
OCnlogn). 

The event point schedule uses 0(n) space and the 
sweep- rectangle status takes 0(1) space. So , the 
total space complexity is 0(n). 

Iasi diM6n3l<ifjnftl sissst 

Now , we extend the above one dimensional approach to 
two dimensions. Here, we have a vertical rectangle of vary- 
ing width d sweeping the plane along the x-axis from -(inf) 
to +(inf). A formal algorithm for this case is given below. 
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Procedure Sueep-rect-tno-dimf 

ilt -finds a closest pair <p,q} of a set of n points in the planel 
begin 

Sort the n points lexicographically first by 
X— coordinate and then by y— coordinate and 
place them in the queue Ef 
S s* Cp $ /» S is the rectangle status *f 
d (inf)} 

If the no, of points in E is more than one do 
begin 

Extract from E the first two elements , find 
the distance betmeen them and assign it to d$ 

S\ - S u ^ > 

athile (E <> <f> ) do 
begin 

Extract from E the first element pi 
Insert pi in S} 

k 

Delete from S all the points having 
abscissae less than (abscissae of pi d) } 
Compute the distance of pi from its nearest 
neighbour candidates in Sf 
Find the minimum of the above distances 
( say dl) } 

Let ql be a point at a distance dl 
from pi} 

d min (dl , d) f 

If d ^ dl then f* update p & q «/ 



begin 


p s- plf 
q s= qlf 
end} 

end} 

end} 

Output d Si (p , q} } 

end} 


Saba gjbra gb wra g 

As In the one dimensional case , here also we can use a 
linear list to implement the event point queue. But the 
data structure for the sweep- rectangle status becomes 
slightly complicated. 

The points in the sweep-rectangle status are kept 
sorted by y-coordinate (the reason will bo clear latter). We 
can use a height balanced tree for storing the points 
.sorted by y-coordinate. An extra link is added to every 
node of the tree. We call this x-link. The x-link of every 
node of the tree is manipulated in such a way that if we 

I 

traverse the tree , using this extra link , we get the same 
points soirtod by their x-coordinates . 

The structure of a node can be Implemented as follows: 
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type 

ptr - ^node; 
node = record 

X t co_ordf 

y s co_prdf 

ylchild 3 ptrf 
yrchild s ptrf 
x—link 3 ptrf 
endf 

Since the event points to be inseirted in the data 
structure comes sorted by their x-coordinates , the insert 
procedure is quite straightforward and can be written as ; 

Procedure IHSERT(p>} 
begin 

1, Insert the nett event point p iP 
the belenced tree , sorted 
by y— coord inetes , 

2- x—link o-f the lest of the 
previously inserted nodes is 
set to point to the currently 
inserted node, 

end} 

LeiiiBa2: The above INSERT procedure takes O(logn) time for 



19 


each point. 

Proof: Stepl of the above INSERT procedure takes O(logn} 

time because insertion in a balanced tree takes 
O(logn) time. Step2 of the procedure takes 0(1) 
time , since it does nothing except set a pointer 
to a node. So, the overall time complexity is 
0( logn) . 

The nodes are traversed through the x-link to execute 
the deletions of the points from the rectangle status. The 
delete operation for the rectangle status ia also straight- 
forward and is as below : 

Procedure DELETE^ 
begin 

t ST ART f 

Pthile t'^.x < (pl.x - d) do 
begin 

delete t#»e nodef 
t t'^. x-link f 

endf 

store t'^.xJink globally to be used as the 

START of delete for next delete operation^ 


endf 
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L«BBa3: The maximum nvimber of dis-bance computations 

required for any newly inserted point is at most 6. 


Proof : 



Let A and B be the vertical sides of the sweep- 
rectangle and p be the newly inserted event point. 
We must find all points q in the sweep-rectangle 
status that are within distance d of p. To do so , 
we note that the most dense packing of points 
inside the sweep-rectangle such that no two points 
are closer than d is as shown in the above figure. 
Thus p has at most 6 nearest neighbour candidates 
which means at most 6 new distances have to be ccaa- 
.puted. 

Theorems: The above algorithm correctly computes a mutually 
closest pair of the planar point set. 

Proof: The proof follows from the observation that at any 

moment of time d and (p,q) gives the distance 
between a closest pair and a closest pair respec- 



•tively of the points seen so far. 

Theore«4: The time complexity of the above algorithm for t-wo 
dimensions is O(nlogn) and the space complexity is 
0(n), 

Proof: To find the nearest neighbour of p we locate the 

node corresponding to p in the sweep-rectangle 
status. This takes O(logn) search time. Then find 
out the successors and predecessors of p (total no. 
of them is bounded above by 6) to get the nearest 
neighbours. This totally takes O(logn) time in the 
worst case. Thus the nearest neighbours of p can 
be found in O(logn) time. Since the total number 
of points in the plane is n , this step takes 
O(nlogn) time. Also insertions into the sweep- 
rectangle status take totally O(nlogn) time and 
deletions takes totally 0(n) time in the worst 
case. The initial sorting stop takes O(nlogn) 
time. 

t- 

Therefore , the overall worst case time complexity 
of the above algorithm is O(nlogn). 

The data structure for the event point schedule 
uses 0(n) space and the sweep-rectangle status also 
uses 0(n) space. So , the overall space complexity 
is 0(n) - 

£. Generalisation higher digeaalonfl 


For simplicity , instead of d-dimensions , let us gen- 
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eralise the algorithm and the data structures to three 
dimensions. Here instead of a rectangle with two vertical 
sides , two hyperplanes parallel to the yz plane , with 
varying distance of separation d , sweeps the three dimen- 
sional space along the x-axis. Initially the points are 
sorted lexicographically first by their x-values , then y- 
values and finally by their z-values , and placed in a queue 
E. The queue E can be implemented in the same way as in the 
one and two dimensional cases above. 

Since the distance between any two points is >= d it is 
seen that here also the maximum number of nearest neighbour 
candidates for every newly Inserted point is 0(1). Hence , 
for each newly inserted point the nvunber of distance compu- 
tations is 6 ( 1 ). The rest of the algorithm and data struc- 
tures are the same as for the two dimensional case. 

# « 

But the optimality of the two dimensional algorithm 

4 

does not carry over to higher dimensions. In higher dimen- 
sions , we don’t find an appropriate data structure for the 
sweep- rectangle status which offers us a O(logn) insertion 
and a O(logn) deletion. Also in the higher dimension the 
simplicity , one of the major advantages of this sweep- 
rectangle technique vanishes.' 

I. aa fC LssioHs 

We have outlined an optimal O(nlogn) algorithm for the 
CLOSEST PAIR problem for "a planar point set , based on the 
novel concept of a sweep-rectangle. However , this optimal- 



ity does not carry over to higher dimensions. But this does 
not detract from the practical usefulness of the algorithm 
since the number of dimensions <= 3 in practical situations. 
It would be interesting to investigate if the sweep-line or 
sweep- rectangle technique can be successfully applied to 
other problems for which optimal algorithms have been 
designed based on other paradigms like the Divide-and- 
Conquer , and vice versa. Guting , for example , has given 
optimal algorithms based on the Divide-and-Conquer strategy 
for various rectangle geometry problems for which optimal 
algorithms based on the sweep- line technique were given ear- 
lier [23 [33. 
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CHAPTER 3 


A STUDY OH THE Kth DEGREE VOROHOI DIAGRAM 


1. I1ITI»D PC ? T I0H 

Given a set S = , qg , ... , q^^} of n points 

(called sites) in the 2-dimensional plane with each point q^ 

represented as an ordered pair (Xj , y^) , 1 = 1 . 2 

n , let d (q^^ , q^) denote the Euclidean distance between 
the two points q^ and q^. The locus of points closer to q^ 
than qj , denoted by h(qj^ , q^) , is one of the half planes 

determined by the bisector BCq^^ , q^) and is = {r | » 

r) < d(qj , r)>. The locus of points closer to q^^ than any 
other point in S , denoted 'by i; (i) , is thus given by V(i) 

= n h(q. , I the intersection of all the half planes 

associated with q^. Vertices of the Voronoi polygon are 
called Voronoi points and their boundary edges are called 
Voronoi edges. The set of Voronoi polygons partitions the 
plane into n regions , some of which may be unbounded , and 
is referred to as the Voronoi diagram V(S) for the set S of 
n points. It has been shown that the Voronoi diagram for a 

set of n points in the plane can be constructed in O(nlogn) 

1 

time. 
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There has been a number of extensioiis and generalise” 

tions of the Voronoi diagram. The definition of the Voronoi 

diagram given above can be easily extended to the L -metric 

P 

where 1 <= p <= (inf) [1] , and the diagram can still be 

constructed in 0(n log n) time , if we allow that the compu- 
tation of the pth root can be done in constant time. 

The second extension consists of considering the Voro- 
noi diagram of a set of objects rather than points. In [23 

the Voronoi diagram for a set of line segments or cii^cles is 

2 

considered and an 0(n log n) time algorithm is given for 
its construction. The time bound was later improved by 
Kirkpatrick to 0(n log n) [3]. 

The third extension focuses on the fact that in the 
Voronoi diagram discussed so far , each polygon is the locus 
of points nearest to one point. To be more precise the 
diagram should be temned the nearest neighbour Voronoi 
diagram. Shamos and Hoey [4] considered the order-k Voro- 
noi diagram of a set of points where each polygon of the 
diagram is associated with k points , k >= 1 , with the pro- 
perty that for any point inside the polygon its k nearest 
neighbours are precisely the associated k points. With the 
order-k diagram the k-nearest neighbours search problem can 
be solved in O(log n ,+ k) time , where the first term 
accounts for point location and the second term for report- 
ing the answer. Properties and a method for the construc- 
tion of the order-k Voronoi diagram can be found in [5] , 
[6] , [4]. 
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The fourth extension is to associate each point with a 

positive weight , resulting in a "weighted" Voronoi diagram 

[7]. The weighted Voronoi diagram consists of n "regions", 

each of which is the locus of points whose weighted distance 

2 

to a given point is minimum. An 0(n ) algorithm for con- 
structing such a weighted diagram can be found in [8]. 

Another extension consists of generalising the diagram 
or triangulation to higher dimensions. Some results in this 
direction can be found in [8] , [10] , [113 , [12]. 

In this chapter we are going to discuss a generalisa- 
tion of the Voronoi diagram , called the higher degree Voro- 
noi diagram. In this higher Voronoi diagram , we partition 
the plane into polvgons in such a way that every polygon 
gives the locus of the points kth closest to a particular 
site. 

Kth degree Voronoi diagram , as it is called , has 
many practical applications. For example , in an informa- 
tion retrieval system , the request of searching for the .kth 
nearest records to a query in a file with n records is quite 
common. And this can be solved using a kth degree Voronoi 
diagram. 

2. DSFIHITIQH QE A Sih 12ESBES TOBQMQJ DIAflBAM 

For a given set S = {q^^ , <l 2 > • • • > ® points in 

the 2“dimensional plane , a kth degree Voronoi diagram , 
denoted by V^ ({q^^ , .... q^}5 , or V^ (S) for short , is a 
partition of the plane into some convex Voronoi polygons. 
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Each kth degree Voronoi polygon , denoted by (<q^ , , 

• . . » <1^> ) is asociated with an ordered set of k points and 

is the locus of the points closest to the site q^^ , second 

closest to the site q.^ kth closest to q^. For a 

particular site we define the corresponding kth degree Voro- 
noi region as 

V (q^j^) = {x I d {x , q^^) is the kth smallest of the 

sequence {d(x , i = 1 , 2 , .. ,n}} 

It is nothing but the union of a disjoint set of Voro- 
noi polygons each of which is associated with an ordered set 
of k points , q^^ being the kth in each. So , V (q^) gives 

the locus of the points kth closest to the site q^^ in the 

plane. V (q^^) may also be empty. 

A kth degree Voronoi diagram can also be expressed in 
terms of higher order Voronoi diagrams. In fact , 

(S) = Vj^ (S) n (S) , 

where Vj^. (S) is the *order-k Voronoi diagram of the 
given set S of n points. 

j3. gfloesgriss qe a Eth degkee moMOi diagram 

In this section , we talk about the general properties 
of a kth degree Voronoi diagram. 

Propl: Each polygon in a kth degree Voronoi diagram is 

associated with an ordered set of k points from the 
set S , and there exists' at most one polygon 
corresponding to each ordered set. 
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Proof: The points inside a polygonal region of the kth 

degree Voronoi diagram is nearest to a point of 
the set S , second nearest to a point qg , third 
nearest to a point q^ and so on. So , we get an 
ordered set (q^^ , qg > ^3 » • • • > of cardinal- 

ity k which identifies the polygon. 

Prop2: In a kth degree Voronoi diagram (S) , if q.is a 

point on a Voronoi edge , then kth nearest neigh- 
bour of q is not unique. 

Proof : Straight forward . 

Props : Total number 4 regions in a kth degree Voronoi 

diagram is. 0{kkl (n-k) ) . 

Proof: The proof comes directly from the fact that the 

total number of regions in the kth order Voronoi 
diagram is 0(k(n-k)). 

A. CQHSTBP.GTIOH Q£ A Kbb CESSES YOS O S O I CIASBAS 

Our algorithm for constructing a kth degree Voronoi 
diagram is very simple and follows the line of Lee [6] for 

the construction of a order-k Voronoi diagram. The algorithm 

12 

is iterative , i.e. , initially we get V (S), then get V 
(S) , and so on until we get V^ (S) from. V^"^ (S) for a 
specified k and a set S of n points. What happens in the 
course of iteration is that the Voronoi polygons of the pre- 
vious degree Voronoi diagram get divided into subpolygons. 
The main difference of the algorithm for the kth order Voro- 
noi diagram from that of the kth degree is that in the 



30 


former case , in "the course of iteration , new points and 


edges 

are 

formed 

but 

in 

the 

process 

some old points and 

edges 

gat deleted. 

But 

in 

the 

latter 

case 

new Voronoi 

points 

and 

edges 

are 

formed 

but no 

point 

or edge gets 


deleted. This fact adds to the complexity of our algorithm 
for constructing a kth degree Voronoi diagram. 

i.l. XBE ALGORITHM 

In this section we shall give an algorithm to partition 
a single Voronoi polygon , say ^ degree 1 

Voronoi diagram V^ (S) to get the subpolygons for V^ (S). 

Suppose that the polygons V^ (<Poi>) V^ (<P .>•) are 

adjacent to polygon V^ {<T?^>) . We want to partition V^ 
(<Pjj^>) into m subregions such that each subregion r^^ is the 
locus of points closer to p^ than to any other points except 

for i = o , 1 , , m-1. To partition V^ thus , 

we have to compute its intersection with V^ ( {p^ , p^^ , ... 

, effect of this is the extension of the edges 

Intersecting at the vertices of V^ interior of 

V^ ( <Pjj^> ) » thereby partitioning the polygon V^ ( <P ^> ) • Let 

. • 1 

the vertices of V (<p^>) be denoted a* » ••• » 

^m-1* Esich edge (1^^ , ^i+1^ polygon is a portion of 

the bisector between p^^ , p^ and is represented by the index 
pair (i , n). By assumption each vertex 1^^ is an intersec- 
tion of three ^dges represented by (i , n) , (i-1 , n) , (i 
, i-1). Lot us denote the edge which is Incident with 

and which is not on the boundary of the polygon V^ (<Pjj>) by 
IN(I^). The following figure shows a typical Voronoi 
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technique. We first obtain the Voronoi diagram for sets of 
three points , P0 » Pjl> J {P2 > P3 » P4> 

extending the edges associated with { IN(Ij^) , } , 

{IN(Ig) , IN(I^)} etc. , respectively , into the interior of 
(<Pjj^>). By merging two adjacent Voronoi diagrams for 
sets of three points , we get the Voronoi diagram for a set 
of six points. Repeating this merge process 1 logg m /,3 } 
times , we will obtain the Voronoi diagram for m points. 
The edges of the diagram which are interior to V^ 
will partition V^ ( <Pj^> ) into m subregions . Fig 2 shows the 
merge process of two Voronoi diagrams for {p^ , P0 , Pj^} and 
{P2 » P3 » 54) 3 shows the final merge process for 

the two Voronoi diagrams for {Pg , P0 , ... » P4} and {pg , 
Pg , ... , Pg} . In Fig 3 , the merge process starts with 
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The technique used to merge two Voronoi diagrams is 
discussed in details in [4] . Here we shall omit the details 
of the B^erge process and describe a method of identifying 
the set of Voronoi points on which the divide “■ and - con- 
quer technique is to be applied. 

Suppose that we have obtained a degree i Voronoi 
diagram V^ (S) , k > i >= 1. We can divide the Voronoi 
points into two groups — ol^l points which exist from previ- 
ous degrees and new points which are just created. It is 



easy 'bo see "bhat only "bhe set of new Voronoi points are 
needed in order to construct the 7^“^^ (S) diagram. Note 
that sill the vertices in (S) are new Voronoi points. Now 
to partition V^ (<Pj^ , Pg , ... , p^>) we first obtain the 
set of new Voronoi points. Such new Voronoi point Ij is 
associated with an edge IN (I.). Now to partition V^ (<p. , 

« j 1 

Pg » ••• » Pj^>) we first obtain the set of new Voronoi 
points and then apply the divide-and-conquer technique to 
the set as just described. After each Voronoi polygon is 
partitioned , we need to associate it with a new set of (i + 
1) points , and at the same time mark the new Voronoi points 
for use in the next iteration. In this manner we can obtain 
the (i + l)th degree diagram V^^^ (S). 

4.2. AH&IiYSIS 

Now let us analyse the running time of the algorithm. 
Suppose that the polygon V^ (<Pj^ ’ ^2 ’ 
partitioned has s new Voronoi points » ^2 ' ' ’ ' ’ ^s‘ 

Since there are O(i'.n) Voronoi polygons in V^ (S) and 0(i n) 
new Voronoi points [63 > the total ntamber of operations 
required is 

^ 0{Sj log s^) = 0(iH log W) 

Since (k-1) iterations are required to obtain a kth degree 
Voronoi diagram , the worst case complexity for the entire 
work is : 
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St. D I MAH I C PFDATIOM QE A Eiii 13EQREE VQBOMOI DIAGRAM 

For on-line applications , dynamic updation of Voronoi 
diagrams is very important. Most of the real-life problems 
demand an on-line algorithm. Since a kth degree Voronoi 
diagram is a very complex data structure , dynamising it 
creates some difficulties. We are not aware of any effi- 
cient dynamising technique for constructing a kth degree 
Voronoi diagram. Here , we propose an easy-to-implement 
dynamising technique. But before we go into the details of 
the algorithm , we’ll like to explore a few more properties 
of a kth degree Voronoi diagram. 

Prop4: Each addition of a point divides an existing 

polygon in at most k parts. 

Proof*. The proof becomes very simple if we think each 
polygon as nothing but an ordered set of points. 
Let a polygon P in a kth degree Voronoi diagram be 
denoted by the ordered set (qj^ , qg , ... i <i.^) , 
Now , the addition of a new point q ^^^^ may divide 
this polygon into sub regions. Since the order of 
(q^ , qg , ... , qj^) remains the same for all of 
the subregions , they’ 11 have these points in the 
same order in their associated ordered sets with 
q^^ inserted between any two consecutive points and 
with qj^ deleted from the set to maintain the cardi- 
nality of the set. -Since there are a maximum of k 
places where q^^^ can be inserted , the property is 
established. 



Prop5: Deletion of an existing point may cause adjacent 

regions to merge. 

Proof : Straight forward . 

Before we come to the formal algorithm , we’ll briefly 
discuss the data stiructure used. 

&.I. DATA STROCTOEE 

We can use any of the standard data structures lilce 
DCEL to implement a kth degree Voronoi diagram. But with 
this Voronoi diagram , we have to store some additional 
information. 

1. Ordered Sets 

With every polygon , we have to store an ordered list 
containing the ordered set associated with that polygon 


2.X>u»ls 

We also have to store the dual of this kth degree 

Voronoi diagram separately which is to be updated with 

the kth degree Voronoi diagram. Initially we shall 

construct a spanning tree of the dual of the Voronoi 

diagram and store it along with the dual. We have to 

update the spanning tree when the dual is updated. 

Construction of the spanning tree from the dual ini- 

2 

tially will require 0(n ) preprocessing time. 


Def : Wo define the dual of a kth degree Voronoi diagram as 
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a graph where each point corresponds to one polygon of 
the kth degree Voronoi diagram and each edge between 
two points gives the adjacency of the corresponding 
polygons. For Example: 



§§S254-.?«Sree_ Voronoi _diagram_For_4_PointSi 
— Voronoi diagram 
— dual 

The dual of a hth degree Voronoi diagram has the fol- 
lowing properties: 

1, The dual may have cycles in it, 

2, The dual may be used to get the neighbours or adjacent 
regions of a given polygon. 

3, A spanning tree of this dual can be used to traverse 

* 

each of the polygons of a kth degree Voronoi diagram 
efficiently mithout repetition . 

Now , we are ready to present the formal algorithms for 
insertion of a new point in the Voronoi diagram and deletion 
of an existing point from it. 



^•2L. imsebtio h 


The following st-eps are taken when a new point 
inserted into the existing Voronoi diagram. 

Procedure Insert <q s new point) s 

Tin *” 

STEPlt Take the dual oY the kth degree 
Voronoi diagram f 

STEP2i Traverse the I'^oronof diagram using 
a spanning tree of it} 

FOR each polygon 00 
BEGIN 

Let the ordered set associated 
with the polygon be (g ^ ^ ^2 r 

... 

So the polygon is represented 
by f g q^)-)f 

Take the bisectors between q 

nn 

and q^ , and q^ , and so 

on and divide the polygon into 
Oik) parts} 

form the ordered set for each 

new subpolygons by inserting 

q in the ordered list (q^ , 

Tin jL 

^2 * ' ^k^ appropriate 

position} 

Oelete- q^ from the ordered lists 
of the subpolygons , if necessary 
to maintain the cardinality of 


9 



3S 


the ordered setf 

update dual and its spanning treep 
END. 


AHALYSIS 

There are 0(k!n) polygons in the kth degree Vomoi 
diagram for a fixed k and S. Each polygon gets divided 
into 0(k) parts. The manipulation of the existing spanning 
tree requires 0(1) time. So the worst case time complexity 
of the above algorithm is 0(k k! n) 

S.a. IKLETIOH 

Deletion of an existing point is not as simple as the 
insertion and creates some major problems. We directly 
give the formal steps for the process of deleting an exist- 
ing point from a kth degree Voronoi diagram. 

Procedure Delete p 

€Let p^ be the point to be deleted! 

STEPl g traverse the Voronoi diagram 
using a spanning tree of the 
dual graph and delete the point 
p. -from all the ordered lists. 

STEP2g Traverse the dual and perform 
the following stepsg 
STEP2.1 g Find neighbours of a 
polygon using the 
adjacency property 



of the dual. 


ST£P2,2s Check if the ordered set 

associated mith any neighbour ing 
polygon is same as that of 
the current polygon , 

If yes , coalesce them 
into a single polygon^ 

STEP2lt There mill be a number of regions 

having an ordered set of cardinal ity 
{k-I}, To' maintain the cardinality 
to k through out the Voronoi diagram 
0 divide those polygons again according 
to the surroundings. Each of the sub- 
polygons generated mill have an ordered 
set of cardinality k, 

CUXVISIOH ALGORITHM HAS ALREADY BEEH 
PRESENTED IN SECTION 41 

STEP4s Update the existing dual and its spanning 
tree! 

AHALT5IS 

Since ■there are 0(k! n) nodes •to'tally , S'bepl 'takes 
0(kl . n) time in the worst case. Let e be the number of 
edges in the dual graph. So the number of Voronoi edges 
also becomes equal to e. Since we traverse each edge of the 
dual graph at most two times , and since number of comparis- 
ons required in Step2.2 for each two ordered sets is 0(k) , 
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Step2 'bakes ■totally 0(ke) time. Since the number of new 
points is 0(k n) , total time required in Step3 is 0{k n log 
n) in the worst case. Dynamic manipulation of the existing 
dual graph and its spanning tree requires 0(n) time in the 
worst case. So the worst case time complexity of the above 
algorithm is 0(max (k! n , ke , k n logn). 

fi. AS mSSSmmhL ALflOSimM EQ& IBE k£h deorks vdboiioi 

Incremental algorithms are very useful for online 
applications. Once we get the insertion algorithm for kth 
degree Voronoi diagram , designing an incremental algorithm 
becomes easy and straight forward. When we start construct- 
ing a kth degree Voronoi diagram incrementally for a fixed k 
, we have to cross through 2 phases. 

PHASEl: 

When the total number of the points currently present 
on the plane < k 

PHASE2: 

When the total number of points present >= k 

Once we reach Phase 2 i.e. the total number points in 
the plane exceeds k , the incremental algorithm becomes same 
as the insertion algorithm discussed in the last section. 
But for Phase 1 , we have to design a separate algorithm. 
For constructing a kth degree Voronoi diagram , we first 
have to cross Phase 1 , then go into the Phase2. While 

inside Phasel we take the following approach: If n be the 



number of points at present on the plane then construct 
(S) i.e. always we'll have a nth degree Voronoi diagram in 
our hand. So the , incremental algorithm , in Phase 1 , takes 
the following form: 

firoctrduref Phase! f 

STEPlf Locate the new point in (S!, 

STEP2t The new point p divides each 
polygon atmost in n parts, 

STEPSs For each polygon , take the bisectors 
betmeen p^ and each of the members 
in its ordered set and use those 
bisectors to divide the polygon in 
0(n> parts. 

STEP4t Update the dual and its spanning tree. 

So after the above steps we get (S’ ) where S’ = S 

0 repeat the above steps till n = k , and then we 

reach Fhase2 , the insertion algorithm. 

The time complexity of the above algorithm in Fhasel is 

tf. 2 

0( Si i!) So the complexity of the entire Incremental 

1-1 - * 

' 2 

algorithm can be calculated to be 0(^i 1! + (n-k)kk! n) = 

0(k^ k! + (n-k)kk! n) = 0(kk! (k + (n-k)n)). 

Z. MSCPSSIQHS 

The algorithm presented here for constructing a kth 
degree Voronoi diagram is iterative and so conceptually 
easy. The strategy proposed here for the dynamic updation of 
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a ktih degree Voronoi diagram is concep'bually easy and also 
easy to implement. The same strategy can also be used for 
the order k Voronoi diagram. There is scope for furthur 
work in the following directions : 

1. Sweep line algorithm for a kth 
degree Voronoi diagram. 

2. Construction of kth degree Voronoi 
diagram in different metrics. 

3. Weighted kth degree Voronoi diagram. 

4. Geodesic kth degree Voronoi diagram 
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CHAPTER 4 


SHORTEST PATH PROBLEM INSIDE A SIMPLE POLYGON IN 

THE PRESENCE OF OBSTACLES 


1. INTBDDP C TJ Q li: 

Recently there has been a significant upsurge of 
results concerning geometry inside a simple polygon which 
include an improved triangulation algorithm [13 » efficient 
algorithms for calculating the geodesic center and the geo- 
desic diameter of a polygon , a number of new shortest-path 
and visibility-related algorithms that require linear amount 
of time beyond a triangulation [2] , and algorithms for link 
distance problems [3]. Some of this work concentrates on 
internal distance analogs of fundamental problems for point 
sets in the Euclidean plane. For example , Toussaint [4] 
developed an algorithm for computing the "relative convex 
hull" of a set of points , which is the shortest cycle con- 
taining all given points and contained in a given simple 
polygon. 

Our present work extends the shortest-path problem 
inside a simple polygon to the case where the simple polygon 
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ontains some obstacles inside it. The problem can be for- 
lally stated as follows: 

Given a 'fixed source poi'nt X i'nside a poly gen P (convex 
or simple} containing a number of obstacles inside it , 
calculate the shortest paths inside P from X to all 
vertices of P and provide a preprocessing of P into a 
^data structure from which the length of the shortest- 
path inside P from X to any desired target point Y can 
be found in time Odog n) ; the path itself can be 
found in time Qdogn k) , nhere k is the number of 
segments along the path. 

Some work has been done previously towards the design 
an algorithm for such preprocessing of P but under the 
sissumptlon that P does not contain any obstacles inside it 
[ 23 . 

Our algorithm uses the familiar sweep-line strategy to 
preprocess P. However , peculiarities of this problem 
necessitate a somewhat non-standard implementation of the 
strategy (it will be evident from the discussions in the 
later sections). 

Possible applications of our algorithm include the 
closest point problem , the nearest post-office problem , 
the problem of finding the shortest-path to any target point 
from a particular source point in the context of a polygonal 
universe , such as an (polygonal) island with interior lakes 
or a polygonal factory floor with interior lawns etc. The 



presence of obstacles inside the polygon makes the problem 
far more useful in real life-applications. 

This chapter is organised as follows. Section 2.1 
describes an algorithm for a convex polygon containing a 
single obstacle in the form of a straight line segment which 
we shall call a line-obstacle hereafter. Section 2.2 
describes an algorithm for the case when P is a simple 
polygon and contains a constant K number of line-obstacles 
inside it , where K >= 1. Section 2.3 describes the algo- 
rithm for preprocessing a convex polygon containing a single 
polygonal obstacle. Lastly in Section 3 we conclude this 
chapter > mentioning the possible extension of the algorithm 
of Section 2.3 to the case when P is simple and number of 
polygonal obstacles inside it is more than one and some 
related open problems. 

Z. AL^RITBM IQ COMEPIS 3HE SHOBTESI EAlfl IEEE QE A SIMEIiE 
mswm HIXB QBSIACtiES IHSIJE 

Let P be a convex or simple polygon with n vertices and 
let s be the given source point on or inside P . For each 
vertex v of P let PA (s,v) denote the Euclidean shortest- 
path f rom. s. to v lying inside P and }PA (s,v)l the length of 
this path. It is well known that ^ PA (s,v) taken over all 
vertices v of P , is a planar tree (rooted at s) , which 

we call the shortest-paih tree of P with respect to s. 
This tree has altogether n nodes , namely the vertices of P 
, and its edges , in this case , are straight line 

segments connecting these nodes to the same point . Our goal 
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is to compute this tree and to partition P in such a way 
that the length of the shortest-path inside P from X to any 
desired target point Y can be found in O(log n) time. 

Before we go to the case of a simple polygon with K 
polygonal obstacles inside , we would like to discuss some 
simpler cases. First we give an algorithm for preprocessing 
P when P is convex and contains a single line-obstacle. 
Then we propose a sweep-line strategy for the second case 
i.e. when P is simple polygon containing K line-obstacles. 
Then wo describe an algorithm when P is convex and contains 
a single polygonal obstacle. Finally we discuss how this 
algorithm can be extended to the case when P is a simple 
polygon and contains more than one polygonal obstacles. 

2.1. COMVEX ESMMl COSTAIHIHG A SlMfiliS LIHE-QBSTACLS 

This case is very simple and we don’t require the 
sweep-line strategy to solve this problem. It is easy to 
see that without any obstacle inside , PA (s,v) is nothing 
but a straight line from s to v. The presence of an obsta- 
cle divides the convex polygon into three sub- regions so 
that shoftest-path to any point in a particular region 
passes through the same corner points. 


hyp- boundary 






In the above figure s is the source point and O is the 
lin©”obstacle inside the polygon P. The straight lines stj^ 
and stg are the tangents to the obstacle 0 drawn from the 
source vertex s. Let the straight lines st^ and stg inter- 
sect the boundary of the convex polygon P at ponts a and b 
respectively. So we get two regions inside the polygon -- I 
, the region [t ^ tg b a] in P invisible from s and V , the 
remaining portion of P visible from s. 

P = V U I 

Now , PA (s,v) is a straight line segment joining s and 
V when v^ V. But PA (s,v) is a polygonal path whose comer 
is the point tj^ or tg when v e I. For v 6 I , 

PA {s,v) = min (stj^ U t^^v , stg VJ tgV). 

The invisible region I is called the obstructed region 
due to the line-obstacle 0 and is denoted by OR (0). We 
can also partition I into regions and Ig so that if v« 

, PA (s,v) = stj U tj^v and if ^*13 ’ = ®^2'^ 

tgV. The locus of the partition boundary is given by : {x : 
jtj^x} - ItgX} = K = Istgi - {stj^i , which is a part of a 
hyperbola and is called the hyp-boundary for the obstacle 0. 
When K > 0 , this hyp-boundary looks like the broken line 

(Fig 1) which partitions I into and Ig. 

The shortest-path to any point in 1^^ from the source 
vertex passes through t^. This point t^ is called the cusp 
of the region 1^ and is denoted by CUSP (I^^). Similarly the 
point tg is the cusp of the region Ig and we can write tg = 
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CUSP (Ig). 

since there are at most 3 sub-regions inside P , any 
point location problem takes 0(1) time. So we can check all 
the vertices of convex polygon and calculate their 
shortest-paths according to the regions where they are 
located in 0(n) time. Also , if an arbitrary target point v 
is given , we can calculate the shortest-path PA (s,v) in 
0(1) time. 

Formally , the above algorithm for partitioning P can 
be written like this: 

Procedure Partition (P s Convax—pol $ s s sour ca- 
ver tax) 

€Let tjj and be the end points of the 
line-obstacle fnsrfde PI 
beifin 

Dr mu tangent lines st ^ and st^ f 

Put a back pointer to s in each point t^ and t^ $ 
Extend st ^ and st^ to meet the boundary of P 
at a and b respectively f 

Divide the region Ct ^ t^ b aJ by the locus 
of V given by /t^vf - ft^vf * K and name the 
si,br,gi<,«s Bbd , 

Let the subregion gives the locus of the points 
mhose shortest-path passes through t 
V P - ( Z^U ; 

CUSP <Z^) f t^ $ 



CUSP ( ; 

CUSP <V} s» 5 ; 

ertUm 


The algorithm for finding the shortest-path of a target 
point from the source point s , after we have obtained the 
shortest-path partitioning of P for s , is given as; 

ProcitUurit f ind-shortesi-^patth CP s Convex-pol $ s z 
Sourcr-point f v g target-point ) 

Begin 

Locate v inside P and find region R 
in ehich v is located $ 

Let X s» CUSP (R) f 
PA (s,v) s* XV U PA ($,x> , nhere 
PA CSfX} is obtained following the 
back pointers from x to s g 
end» 


AMALTSIS 

It is very easy to see that above two algorithms take 
0(1) time each. 

Z.Z. SIMPLS m TOOB OQHTMHIHS K IiIHE QBSTACbES 

This problem is far more complicated than the previous 
one. So before we go to the main problem we'll explore the 
case when polygon P Is convex. We use the familiar sweep- 
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; i. T.. KANPUR 


icc. No. a.10b4.1.3.9 

li.no 'tochnlQUo lior© iio portition "fc!!© polygon P. ThOr© © 3 r© © 
cons'tan'b K nuiaber of line-obstacl©© Inside "the polygon. W© 
can define a "type of ordering of these obstacles Inside the 
polygon. 

Def: We assign a level number to each of the obstacles. 
The level of an obstacle can be defined as follows. We 
define the level of an end point v of an obstacle as 
one more than the level of the last corner point of the 
ahortest-path from the source point s to v. The level 
of the source point is taken to be zero. The level of 
a line obstacle is equal to the maximum of the levels 
of its end 3 E>oints. So , level of an obstacle visible 
from the source point s is equal to one. 

We can use a sweep- line technique to find out the level 
of an obstacle. In this technique , a line sweeps from the 
source vertex towards the obstacles and stops outside the 
polygon P. In the course of the sweep it calculates the 
levels of the obstacles it meets inside the polygon. Con- 
currently > it partitions the polygon P into shortest-path 
regions. .Instead of the sweep-line status , here we use 
PARTITIONIST ATOS.‘ At any moment of time the PARTITION_STATUS 
gives the partitions of the polygon P created till that 
time. As shown in the last section , for each obstacle 0^ , 
we get an obstructed region OR (0^) . In our algorithm , 
when the sweep-line meets an line-obstacle , its obstructed 
region is found out by constructing the tangent lines to its 
end points. Every end point has two fields as given by the 



type declaration below. 


Type 

pointer » '^en<l_poxnt 

end__point ® record 

path f pointer 

length s integer 

end $ 

When the sweep- line meets an end point , it sets the 
pointer field of the end point to point to the origin of the 
tangent line to this point. Also it finds out the 
shortest-path length to the origin of the tangent line , 
calculates the 8hortest_path length to the end point and 
stores it in the specified field. 

For constructing a tangent line to an end point , the 
algoribhu locates the point in the present partitioning as 
given by the PAKTITION.STATOS and finds out the region it is 
located. The CUSP of the above region is taken to be the 
origin of the tangent line to the end point. 

After the obstructed region OR (0^^) for a line-obstacle 
0^ is found out I we partition the region in two parts by 
the corresponding hyp-boundary , as discussed in the last 
section. OR (0^^) is then inserted in the PARTITIOH.STATUS. 
Some intersections may be possible between more than one 
obstructed regions .in the PARTITION.STATUS. So . 
PARTITION_STATUS is to be updated accordingly. Before we 
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present the formal algorithm , we’ 11 give a list of observa 
tions f required for the updation of the PARTITION_STATUS : 

Obsl: A hyp-boundary vanishes when it comes out of th 

corresponding obstructed region. 



0bs2: Hhen tmo obstructed regions intersect , their hyp— 

bounder ies may also intersect nith each other. Inter- 
section of tmo obstructed regions may require construc- 
tion of ne» hyp— bound ar ies as shown in the following 

figuret 



ObsS-. X « * /.y(.-l.«und.ry 



if xt Meets a line-obstAcle as shorn in the follottinq 
figures 



0l>s4: A tamgent line terminates mhen it meets a tangent 

line of an enclosing obstructed region. An obstructed 
region is called to be enclosing for a giVen obstacle 
0^ if the shortest path to the obstructed region of 0^ 
passes through an end point of the enclosing obstacle 
(shomn in the following figure!. 
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Formally our algorithm can be given as follows; 

Procedure Wafre-partition (P t ConvBx_pojlJ> f 
Initial isatioR s 
Begin 

CliECKE1>_STACK f f 

PARTIT J9H_JST ATiiS i The polygon boundary edges Mhich 

intersect a vertical line through 
the source vertex J f 

Sort the vertices of P and the end points 
of the linm--obstacles along the abscissa 
and store them in a queue E , the event 
queue for the sweep-line f 
end f 
Segi n 

Sweep a line from the source vertex s to v-daf.) 
along the abscissa t 

For each event point e^ perform the following stepsf 
STEPlt If e^ is one of the end points 

of a line-obstacle , take the following 
actions t 

STEPlmls Locate e^ in the PARTITION_STATiiS ; 

Let the region e^ is located be R. $ 

STEPl,2s p^ Ci/SP (R\) $ 

ST£P1.3$ Join points p^ and e^ and extend 
the lipe, 

let the line be denoted by TL CeJ, 

Me call this line as tangent line for the po 



♦ Insert the line in the 

Pi^TJTIOH_STATUSmCheck if the line inter sects nit I 

mny of the polygon boundary edges in the 

P ART IT I OH ^ST AT US m Find out the inter sect ions » 

STEPl m4s If 1 be the level of the line-obstacle nhose 

one of the end points is p^ , then 

assign level 1 -f- 1 to the point e. s 

‘1 

STEPl,5s If e^ be an end point of a line— obstacle 
mhose other end point Oj has already been 
processed and if p. be the origin of 
the tangent line TL <ejl , then find 
the locus ij. of the point x given by 
the equations 

/PA is , e .>/ * /e . xl /PA (s , e .) / * /e. xi 
This is «i» equation of a hyperbola and ne 
call this curve 1 . as hyp-boundary nith 
focii .»<( 

Extend 1^ to meet the bondary of 
P or TA (e^) or TA <e after nhich 
it vanishes • TA (e^ ) , TA 

form imo regions. Let the region 
bounded by TA and 1^ be denoted 

by I^ and the other region bounded by 
TA denoted by 

Then 0 

CUSP (1^1 s*e^ 


CUSP <I.> i* e^ 



Insert ihw mtutve three regions in the 
P^VtTITION^STATUSM Check intersection nith 
the other regions present in PARTIRION__STATUS 
end eccordingiy modity PARTITIOH_STATUS.This 
modlificetions should he done in line of the 
observations listed above. 

The level of the obstacle z- eaximum of the 
levels of its end points. 

STEP2t If the event point e^ is one of the vertices 

of the polygon P , then take the follotzing actions z 
STEP2.is Remove the boundary edge of P mhich ends 
at the point e^ from the PARTITIOH_ 

STATUS. Remove mith this edge all the lines 
intersecting it from the PARTIT10H_STATUS f 
STEP2.2S Insert the nem boundary edge of P mhich start 
from the event point e^ in the PARTITIO_ 
STATUS. Check if any line , currently inside 
PARTITIOH^STATUS , intersects this ne» edge, 
find out all the intersections f 


Onoo th« above algoritha ia known it is not a major 
probleii to extendi it to the case when P is a simple polygon 
with a veirtioes. The only problem in the case of a simple 
polygon is that the polygon boundary also sometimes behave 
as an obstacle. As seen in the following figure , the obs- 
tacle 0 is not visible from the source vertex. So con- 
structing a tangent line poses some difficulties. We have 



■bo take a polygonal path to reach the end points of tlie 
line-obstacle 0. 



We can circumvent this difficulty by performing a bit 
of preprocessing. What we do is that before we apply our 
sweep-line algorithm , we perform a shortest-path partition- 
ing of the simple polygon P , removing all the obstacles 
from inside. There is already an 0(n log log n) algorithm 
available [2} to do this type of partitioning of a simple 
polygon without any obstacles inside. We can use this algo- 
rithm as a preprocessing algorithm. Once we do this initial 
partitioning , we can locate any point in P in 0(log n) time 
, find the shortest path , draw the tangent lines and modify 
the pairbitioning by the sweep-line strategy. 

Once the partitioning is available , finding a shortest 
path from the given source point to an arbitrary target 
point becomes an easy job. The same f ind_sho3rtest_jpath 
algorithm , presented in the last section , can also be used 
here for that purpose. 


Before we start analysing the time complexity of the 
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above algorithm we shall like to mention that for some dis- 
tributions of obstacles and the source point it may be 
required that the above algorithm be applied twice. Once 
while sweeping the polygon from the source vertex to +(inf) 
and nekt time sweeping from the source vertex to -(inf). An 
example of that type of distribution is given below: 



AHAItYSIS 

It is easy to see that the number of partitions of th< 
polygon P is totally dependent on the constant K and th< 
total number of levels of line-obstacles , which is als< 
constant for a distribution of line-obstacles and a givei 
source point. So , the time complexity of the above sweep- 
line algorithm is 0(1) when P is convex. But when P is i 
simple polygon , we perform an initial preprocessing whic] 
partitions P in 0(n) regions. For each end point of th 
line-obstacles we have to do a point location. Since ther 
are 0(1) line-obstacles inside the polygon , the sweep-lin 
algorithm takes 0(log n) time in the worst case. For conve 



polygon , the f ind_shortest_path algorithm takes 0(1) time 
and for a simple polygon it takes 0(log n) time. 

2-3. QfMl E X. EQ Ii YQOM CX?HTAIH I H S A SJHgliE FOLYBCMat. oBSTAnr.TC 

t 

A typical partitioning of the convex polygon P is given 
below when there is a single polygonal (convex) obstacle 
inside it. 



The algorithm is easy and straight forward in this 
case. We don’t require any sweep-line strategy to solve 
this problem. The Steps required to do this partitioning is 
given below: 


STEPl; Draw two tangents from the source to the 
polygon P. Let the points of contact be 
t^ and t^. Extend the 
tangent lines to meet the polygon 
boundary at a and b respectively. 

The region £t > ^2 ^ 

is invisible from s and we denote 


this region by I 


inside 


STEPS: Extend the 
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the region I. Let the region formed be 

«> 

called as For any point in 

the region , the shortest-path 

passes through either or t^ ; 

STEPS’* Divide each region in two parts 
by the hyp-boundary 1^^ with focii 
tj.i and (as dlsouasad 

earlier). The hyp-boundary 1^^ 
divides the region into regions 

®il **12- 

STEP4; Lot the shortest-paths for the points in 
Rj^j and P®®® through 

tj^j and t^ respectively. 

Then , 

CDSP (R^i) == 

(aJSP (Rj^g) != \ 

Once we perform the partitioning finding a shortest- 
path from s to an arbitrary target point X is just a point 
location problem. The same find_shortest_path algorithm can 
be used here* for that purpose. 

Once we have an algorithm for a convex polygon we can 
extend it to the case when P is a simple polygon in the same 
fashion described in the last section, 

2. PI8CT8SX0BS. 

We can extend the above algoritim for a single polygo- 
nal obstacle to the case when the hhmber of polygonal obsta- 



cles is more than one. We can use the same sweep-line tech- 
nique for this purpose. Since in this case the partitioning 
of the polygon becomes very complicated , no formal algo- 
rithm is given here. 

This problem can also be solved by taking each polygo- 
nal obstacle as a combination of line-obstacles and then 
using the same sweep-line algorithm. 

It remains an open problem to deteirmine whether such a 
partitioning is possible for a three dimensional case. Also 
there are a number of open problems closely related to this 
shortest -path problem. A problem of major importance is to 
see whether it is possible to design an algorithm for geo- 
desic Voronoi diagram when the simple polygon contains a 
number of obstacles. 
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CHAPTER 5 


CONCLOSIOH 


In l/his i/hosis wo havo looked ali various kinds of 
proximity problems , 

In chapter 2 we have outlined an optinal 0(n log n) 
algorithm for the closnsi pmir problem for a planar point 
set , based on the novel concept of a sweep-rectangle. How- 
ever , this optimality does not carry over to higher dimen- 
sions. But this does not detract from the practical useful- 
ness of the algorithm since the number of dimensions is less 
than three in most practical applications. The sweep-line 
strategy is primarily suited for problems regarding lino 
segments (for example : line intersection problem). The 
strategy called sweep- rectangle technique , proposed in this 
chapter , is a variation of sweep- line and is suitable for 
problems regarding point sets. It would be interesting to 
investigate if the sweep-rectangle technique can bo success- 
fully applied to other problems regarding points. 

In chapter 3 we have discussed a type of generalisation 
of Voronoi diagram called nth degrmm Voronoi diagram. In 
this chapter we have also proposed a new algorithm for 
dynamic updates of a Kth degree Voronoi diagram. It remains 
an open problem to determine idiethar a specialised dynamic 



data structure can be designed for the Kth degree Voronoi 
diagram. There is scope of further work in the following 
directions: dynamisation of order-K Voronoi diagram ; dynam- 
isation of geodesic Voronoi diagram ; designing an algorithm 
to construct order-K or Kth degree Voronoi diagram ; design- 
ing sweep-line algorithm for geodesic Voronoi diagram ; Kth 
degree Voronoi diagram , and order-K Voronoi diagram. 

In the third chapter we have dealt with a different 
kind of proximity problem called shortest-path partitioning 
of a simple polygon. In this chapter we have proposed a 
sweep-line algorithm to partition a simple polygon contain- 
ing a number of obstacles. It remains an open problem to 
determine whether such a partitioning is possible for three 
dimensional case. An open problem of major importance is to 
see whether it is possible to design an algorithm for geo- 
desjfc Voronoi dimgrmm when the simple polygon contains a 



